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We study the scattering of small color dipoles (e.g., heavy quarkonium states) 
at low energies. We find that even though the couplings of color dipoles to the 
gluon field can be described in perturbation theory, at large distances the interac- 
tion becomes totally non-perturbative. The structure of the scattering amplitude, 
however, is fixed by the (broken) chiral and scale symmetries of QCD; the lead- 
ing long— distance contribution arises from the correlated two-pion exchange. We 
use the spectral representation technique to evaluate both perturbative and non- 
perturbative contributions to the scattering amplitude. Our main result is the sum 
rule which relates the overall strength of the non-pcrturbative interaction between 
color dipoles to the energy density of QCD vacuum. 



1 Introduction 

In a 1972 article entitled "Zero pion mass limit in interaction at very high 
energies" 0, AA. Anselm and V.N. Gribov posed an interesting question: what 
is the total cross section of hadron scattering in the chiral limit of — > 0? 
On one hand, as everyone believes since the pioneering work of H. Yukawa, the 
range of strong interactions is determined by the mass of the lightest meson, 
i.e. is proportional to ~ rn^ 1 . The total cross sections may then be expected 
to scale as ~ m~ 2 , and would tend to infinity as 0. On the other hand, 

soft-pion theorems, which proved to be very useful in understanding low- 
energy hadronic phenomena, state that hadronic amplitudes do not possess 
singularities in the limit m„ — > 0, and one expects that the theory must remain 
self-consistent in the limit of the vanishing pion mass. 



a Invited talk at the Third "Continuous Advances in QCD" Workshop dedicated to the 
memory of V.N. Gribov; Minneapolis, April 16—19, 1998. 
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At first glance, the advent of QCD has not made this problem any eas- 
ier; on the contrary, the presence of massless gluons in the theory apparently 
introduces another long-range interaction. In this article, we try to address 
this problem considering the scattering of small color dipoles. The reason for 
choosing this, somewhat specific, example is simple: asymptotic freedom dic- 
tates that strong coupling becomes weak at short distances, and since the small 
size of dipoles r introduces a natural infrared cut-off r <C Aq^ d , one expects 
that their interactions can be systematically treated in QCD. This approach 
was pursued both at lowiSBIu and high@ energies. 

It has been found, however, that at sufficiently high energies the pertur- 
bative description of "onium-onium" scattering breaks downifl The physical 
reason for this is easy to understand: the higher the energy, the larger impact 
parameters contribute to the scattering, and at large transverse distances the 
perturbation theory inevitably fails, since the virtualities of partons in the lad- 
der diffuse to small values ("Gribov diffusion"). At this point, the following 
questions arise: Does this mean that the problem becomes untreatable? Does 
the same difficulty appear at large distances in low-energy scattering? And, 
finally, what (if any) is the role played by pions? 

2 Perturbation Theory 

The Wilson operator product expansion allows one to write down the scat- 
tering amplitude (in the Born approximation) of two small color dipoles in the 
following fornH: 



where x = (t, R), Oj(x) is the set of local gauge- invariant operators expressible 
in terms of gluon fields, and are the coefficients which reflect the structure 
of the color dipole. At small (compared to the binding energy of the dipole) 
energies, the leading operator in ([!]) is the square of the chromo-electric field 
(l/2)g 2 E 2 !{] — other twist-two operators contain covariant derivatives leading 
to the powers of momentum in the amplitude and are therefore suppressed at 
small energies. 




(1) 
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Keeping only this leading operator, we can rewrite (|IJ) in a simple form 
V(R) = -i(d 2 ^) 2 J dt(0\T ~g 2 E 2 (0) ±g 2 E 2 (t,R)\0), (2) 
where d-i is the corresponding Wilson coefficient defined by 

e 3iV H a + e 

where we have explicitly factored out the dependence on the quarkonium Bohr 
radius a® and the Rydberg energy eo; N is the number of colors, and \<p) is 
the quarkonium wave function, which is Coulomb in the heavy quark limit^ 
The factors eto and r ~ 1/eg represent the characteristic dimension and fluctu- 
ation time of the color dipole, respectively. The approximate expression (||) is 
justified when the wavelength of gluon fields is large compared to ao and they 
change slowly compared to 1/eo. 

In physical terms, the structure of (||) is transparent: it describes the 
elastic scattering of two dipoles which act on each other by chromo-electric 
dipole fields; color neutrality permits only the square of dipole interaction in 
the elastic scattering. 

The amplitude (|^) was evaluated before! in perturbative QCD using func- 
tional methods. For our purposes, however, it is convenient to use a different 
technique based on the spectral representation. As a first application of this 
approach, we will reproduce the result ofi by a different method. 

It is convenient, as a first step, to express g 2 E 2 in terms of the gluon field 
strength tensor El: 

ff 2 E 2 = -\g 2 G^G^ + .g 2 (-G 0Q G« + \g m G a ^) = + g 2 9^ ] (4) 



4 

where 



Note that as a consequence of scale anomaly 0, 8^ is the trace of the energy- 
momentum tensor of QCD in the chiral limit of vanishing light quark masses, 

6 The Wilson coefficients J21 evaluated in the large N limit, are available for s| and pQ 
quarkonium states. 



3 



and the (3 function in (||) docs not contain the contribution of heavy quarksQ 
{i.e. &=|(11JV-2JV» = 9). 

Let us now write down the spectral representation for the correlator of the 
trace of energy-momentum tensor: 

(O|T0£(O)0£(z)|O> = / da 2 p e (<J 2 )A F (x;a 2 ), (6) 
where the spectral density is defined by 

p e (k 2 )0(k o ) = J2(^) 3 S\ Pn - fc)|(O|0»| 2 , (7) 

n 

and 

iA F (x;a 2 ) = ij^6(k 2 -a 2 )e(k a )(e- lkx 6(x°) + e lkx 8(~x )) (8) 

is the Feynman propagator of a scalar field. 
Using the representation @ in (^), we get 

v,(R) = -i(dA)\l^-) 2 j dt j dSe,{a 2 )&Fix;c*) 

We see that the potential (||) can be represented as a superposition of Yukawa 
potentials corresponding to the exchange of scalar quanta of mass a. 

Our analysis so far has been completely general; the dynamics enters 
through the spectral density (Q). Let us first evaluate this quantity in pertur- 
bation theory. In this case we have to evaluate the contribution of two-gluon 
states (see Fig. 1(a)) defined by 

pf{q 2 ) = ^(27r) 3 5 4 (p 1+J32 - (? )|(0|Cbieia,P2e 2 6)| 2 , (10) 

where the phase-space integral should be understood as well as the summations 
over the polarization and the color indices of the gluons. 

The evaluation of (|l0|) is straightforward: for SU(N), one has 

pt 2 ( bg 2 \ 2 N 2 -1 4 
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Figure 1: Contributions to the scattering amplitude from (a) two gluon exchange and (b) 
correlated two pion exchange. 



Then, substituting ( pi] ) into (^) and performing the integration over invariant 
mass a 2 , we get, for N = 3 



WJ0 . V(i5) 2 ^. (12) 

The oc R~ 7 dependence of the potential ([l2|) is a classical result known 
from atomic physics 111; as is apparent in our derivation (note the time inte- 
gration in (||)), the extra Rr 1 as compared to the Van der Waals potential 
oc R~ 6 is the consequence of the fact that the dipoles we consider fluctuate in 
time, and the characteristic fluctuation time r ~ 6q , is small compared to the 
spatial separation of the "onia" : r <C R. 

Finally we note that in perturbation theory 9fi is of order g 2 , (even though 
the matrix element of the 6fi, as will be discussed below, is in general non- 



perturbative), and accordingly the potential (12) has the prefactor g A . Then 
the second term in (Q) gives the contribution of the same order in g; this 
contribution is due to the tensor 2 ++ state of two gluons and can be evaluated 
in a completely analogous way. Adding this contribution to (|l2|) , we reproduce 
the result of Ref.i: 

v ^ = <^t)^h < 13 » 

note that our d 2 is related to the d-2 of Ref. i by ^doeo = dig 2 . 
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3 Beyond the Perturbation Theory: The Role of Pions 

At large distances, the perturbative description breaks down, because, as can 
be clearly seen from (^), the potential becomes determined by the spectral 
density at small q 2 , where the transverse momenta of the gluons become small. 
To see this explicitly in the dispersive language, let us consider the correlator 

U(q 2 ) = J dW«*<O|T0£(zWO)|O> = J do 2 J Q } U P_^ (14) 

An important low-energy theorem@ states that, as a consequence of the broken 
scale invariance, 

n(0) = -4(0|^G Q/3a G° s (0)|0). (15) 
2# 



The operator on the r.h.s. of ([T^) is regularized by subtracting the contribu- 
tion of perturbation theory. The vacuum expectation value of this operator 
therefore measures the energy density of non-perturbative fluctuations in the 
QCD vacuumlll. The low-energy theorem ( |l5|) therefore implies the sum rule 
for the spectral density El: 

/ ^[Pt Y V 2 ) - pgV)] = -4(01^0^^10) = -16 £vac + 0, (16) 

where p^ t {a 2 ) is given by (O), and the vacuum energy density is e vac = 
(1/4) (0£) ~ -(0.24 GeV) 4 ElTsince the physical spectral density, pP hys , should 
approach the perturbative one at high a 2 , the integral in ([!(]) can accumulate 
its value required by the r.h.s. only in the region of relatively small a 2 ^. 

At small invariant masses, we have to saturate the physical spectral density 
in the sum rule (|l6| ) by the lightest state allowed in the scalar channel - two 
pions: 

pf y V) - ^(27r) 3 <5 4 (p 1 +p 2 -q)|(0|^k(pi) 7 r(p 2 ))| 2 , (17) 

where, as in (JlC|) , the phase-space integral is understood. 

Since, according to (g), 6% is gluonic operator, the evaluation of dl7j ) re- 
quires the knowledge of the coupling of gluons to pions. This is a purely 

c The analysis of sum rules shows_however that the approach to the asymptotic freedom in 
the scalar channel is rather slowLD. 
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non-perturbative problem, but it can nevertheless be rigorously solved, as it 
was shown in Ref.S (see alsoEl). The idea of 111 is the following: at small pion 
momenta, the energy-momentum tensor can be accurately computed using the 
low-energy chiral Lagrangian: 

0* = -d^dV + 2mlTT a TT a + ■■■ (18) 

Substituting this expression into (|l7|), in the chiral limit of vanishing pion mass 
one gets an elegant result @ 

{Q M G "f>*G% fl \n + *-) = <?. (19) 

The result (|l9| ) can actually be generalized for the coupling of gluons to 
any number of pions. Indeed, consider the Lagrangian of non-linear a model, 

f 2 

C = i*- tr d^Ud^ - A tr (MU + UH^) , (20) 

where U = exp (2iw/f w ), ir = ir a T a , tr T a T b = \5 ab , M is the quark mass 
matrix, and A = m^f^/m, with rh being the average light quark mass. Using 
the mathematical identity for a generic matrix A, 



^[exp(A)] = f ds exp(sA) d^A cxp((l - s)A), 
Jo 



(21) 



one can explicitly evaluate the trace of the energy-momentum tensor for the 
Lagrangian (|20|), with the result 

f 2 

6» = -2 J -f tr d^Ud^ + 4A tr (MU + U^M^) . (22) 



Expanding (|22j) in powers of pion field, one can generalize ( |18| ) for any (even) 
number of pions; to lowest order, we reproduce (|l8|). 

Now that we know the coupling of gluons to the two pion state, the pion- 
pair contribution to the spectral density (17) can be easily computed by per- 
forming the simple phase space integration, with the result 

<FM = nb*- (23) 
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Multi-pion contribution can be evaluated using (|22|); its influence will be dis- 
cussed elsewhere. However the dominant contribution at small invariant masses 
a, in which we are primarily interested here, comes from the tttt state. 

Coming back to the initial expressions (||) , (|) we find that to complete our 
derivation of the scattering amplitude we need to evaluate also the transition 
matrix elementEl of the second term in (|J), (0\g 2 9^ \irir). This tensor operator 
was discussed in the previous Section, where we have found that it contributes 
a substantial fraction, 8/23, to the complete perturbative result. However, 
unlike the scalar operator, the tensor term is not coupled to the anomaly, and 
therefore {0\g 2 9^ \tttt) ~ g 2 , where the coupling has to be evaluated at the 
heavy quarkonium scale. This contribution therefore is of higher order in the 
strong coupling, and will vanish in the heavy quark limit. Omitting it, we 
come to the following low-energy expression, 

q2, \ _ ( 8tt 2 \ 2 . r,, 2 



(fVE^M = ^ — 1 q 2 + 0{a s ,mi). (24) 

Thus, in the heavy quark limit, the matrix element in question is known up to 
a s and m\ corrections. 

The result (|24| ) has been derived in the chiral limit; the most important 
correction coming from the finite mass of the pion is the phase space threshold; 
we correct for it by writing down the spectral density in the form (q 2 > 4m 2 ), 

which should be valid at small q 2 . Substituting ( p5| ) in @ , we get the potential 
due to the 7T7t exchange; at large R 



The same functional dependence of tttt exchange at large R has been obtained 
previously in Ref.i, but up to an unknown constant; in our approach, the over- 
all strength of the interaction is fixed. Note that, unlike the perturbative result 
(O) which is manifestly ~ g 4 , the amplitude (Eq) is ~ g° - this "anomalously" 



strong interaction is the consequence of scale anomaly^ 

The low-energy theoremllllll not only allow us to evaluate explicitly the 
contribution of uncorrelated tttt exchange; they also tell us that this contribu- 
tion alone is not the full story yet. Indeed, the tttt spectral density (25) alone 
cannot saturate the sum rule (|l7J) - at high er 2 , the physical spectral density 
approaches the spectral density of perturbation theory, so the integral in (16) 
does not get any contribution; at small a 2 , the tttt spectral density (p5|), ac- 
cording to the chiral and scale symmetries is suppressed by ~ cr 4 . The low 
energy theorems require the presence of resonant enhancement (s) El in the ++ 
7T7T, and perhaps multi-pion, KK and r\r\ channels as well. Here we will leave 
this interesting problem aside, and study only the influence of these resonances 
on the potential between the color dipolcs. To do this, we introduce the pion 
scalar form factor F(q 2 ) and write down the spectral density as 

q 2 -4r?4\ 1/2 



tfV) = ^{^^) (27) 

The form factor is directly related to the experimental tttt phase shifts by the 
Omnes-Muskhelishvili equation with the solution 



F(t) = exp 



t[ ds - 



(28) 



s(s — t — ie) 

with this formula we can make a full use of the experimental information on 
the tttt correlation. For our calculation we have used a simple analytic form@ 
for the phase shift 5® which was shown to fit the experimental data up to 
s W7r ~ (1.2 GeV) 2 . There are two main contributions to the spectral density, 
pg v , which may be interpreted as the broad a and narrow fo resonances. 

In Fig. U we show the resulting potential between two J/ip states. In 
computing it, we assumed that Coulomb relations for the Bohr radius and the 
Rydberg energy ao = 4/(3 a s m) and eo = (4/3 a s )~ 2 m = 1/a 2 , m hold for the 
J/ip. Using as an input e = 2Afo - M(J/^)=642 MeV and m=1.5 GeV, we 
get a s =0.87 and a =0.20 fm. 

It can be clearly seen from Fig. |^ that at large distances the non-perturbative 
interaction dominates over the perturbative one. To evaluate the amplitude, 

d Ol course, in the heavy quark limit the amplitude ( prj| ) will nevertheless vanish, since ao — » 
and eo — > oo. 
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Figure 2: The potential between two J ftp's (bold solid line); the perturbative contribution 
V pt (dashed line) was evaluated within the invariant mass range <r > 2 GeV in the spectral 
density; V np (long-dashed line) is the non-perturbative contribution. 

we had to use an experimental input - the tttt phase shifts, and the detailed 
shape of the potential depends on this input. However, as we will now see, the 
dominance of the non-perturbative interaction is a model-independent con- 
sequence of the low-energy theorems. Moreover, its overall strength can be 
shown to be completely determined by the energy density of non-perturbative 
vacuum of QCD. 

4 The Sum Rule 

Consider the integral over the non-perturbative part of the potential; according 
to (^), it can be written down as 

/>oo 

/ d 3 R(V{R) -V pt {R)) 

J a 
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= -(d 2 ^) 2 (^) 2 /^(^ 2 )-p pt (- 2 ))r(2,-)- (29) 

As a consequence of asymptotic freedom, p(a 2 ) — > p pt (cr 2 ) at high cr, so the 
integral in ( p9| ) attains its value in the region of a < <Jo, where <7 is some 
characteristic scale at which the perturbative regime sets in. In the heavy 
quark limit the size of quarkonium a ~ {a s m) — > 0, and when eroa *C 1, 
r(2, era) ~ 1 in the entire region of integration. The integral in then, 
up to perturbative corrections ~ 0(g 4 ), coincides with the integral in ([l6|). 
Therefore we can re-write (p9|), in the heavy quark limit, as a sum rule 

// _ a 2 \ 2 /Air 2 \ 2 
d 3 R(y(i?)-^ 4 (i?)) --(d 2 J) (-£-) 16|e„«c|, (30) 

which expresses the overall strength of the interaction between two color dipoles 
in terms of the energy density of the non-perturbative QCD vacuum. 



5 Final Remarks 

What are the implications of our results? First, the pion clouds which dominate 
interactions of small color dipoles at low energies, as revealed by our analysis, 
may as well be important in high energy scattering; this was suggested long 
time ago on general grounds BEJ. However it is not yet clear if one can extend 
our approach to the scattering at high energies^]. 

Second, the fact that pions (and therefore light quarks) dominate the long- 
distance interactions of heavy quark systems is important for the lattice QCD 
simulations. Our findings suggest that to extract the information on the prop- 
erties and interactions of heavy quarkonia from lattice QCD one should go 
beyond the "quenched" approximation. 
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